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ABSTRACT 



The kinematical evolution of four EUV waves, well observed by the Extreme 
Ultraviolet Imager (EUVI) onboard the Solar- Terrestrial Relations Observatory 
(STEREO), is studied by visually tracking the wave fronts as well as by a semi- 
automatized perturbation profile method leading to results matching each other 
within the error limits. The derived mean velocities of the events under study 
lie in the range of 220-350 kms~^. The fastest of the events (May 19, 2007) 
reveals a significant deceleration of ~ —190 m s~^ while the others are consistent 
with a constant velocity during the wave propagation. The evolution of the 
maximum intensity values reveals initial intensification by 20 up to 70%, and 
decays to original levels within 40-60 min, while the width at half maximum 
and full maximum of the perturbation profiles are broadening by a factor of 
2-4. The integral below the perturbation profile remains basically constant in 
two cases, while it shows a decrease by a factor of 3-4 in the other two cases. 
From the peak perturbation amplitudes we estimate the corresponding magneto- 
sonic Mach numbers M^^ which are in the range of 1.08-1.21. The perturbation 
profiles reveal three distinct features behind the propagating wave fronts: coronal 
dimmings, stationary brightenings and rarefaction regions. All of them appear 
after the wave passage and are only slowly fading away. Our findings indicate 
that the events under study are weak shock fast-mode MHD waves initiated by 
the CME lateral expansion. 

Subject headings: EUVI waves - perturbation profiles - corona - stationary brightenings 
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Introduction 



Large-scale, larg:e-ampli tude disturb a nces 



have been first observed by 



Moses et al. 



propa gating through t h e sola r corona 



fll997h and 



Thompson et al.l (11998) in images 



recorded by the Extreme Ultraviolet Imaging Telescope (EIT; 



onboard the Solar and Heliospheric Observatory (SoHO; 



Delaboudiniere et al 



Domingo et al. 



19951 ) 



19951 ). thereafter 



called 'EIT waves'. They were originally interpreted as t he coronal count erparts o 



chromospheric Moreton waves (iMoreton fc Ramsey 



19601 ) as suggested by lUchidal (119681 ) 



the 



with his coronal fast-mode MHD wave interpretation of Moret on waves. Though this 



Warmuth et al. 



interpretation has been confirmed for severa. 



Vrsnak et al. 



2001 



2006 



Pohiolainen et al. 



Veronig et al. 



2006 



2001 



case studies 



Vrsnak et al. 



Muhr et al. 



e.g. 



Thompson et al. 



2002 



1999 



Khan fc Aurass 



2002; 



20101 ). the statistical characteristics of 



EUV and Moreton waves are quite different. In particular, EIT waves generally propagate 
at much lower velo c ities, v ~ 200-400 kms~^, than Moreton waves, v ~ 1000 kms~^ (e.g. 



Klassen et al. 



2000 



Biesecker et al 



2002; 



Thompson fc Myers 



20091 ) . Thus, there is still an 



ongoing debate of whether EIT waves are really the coronal counterparts of Moreton waves 
Additionally, it is still an open discussion whether they are caused by the flare explosive 



energy release or t 



2004 



Oliver et al. 



he erupting OME (e.g . 



2005 



Vrsnak fc Oliver 



Warmuth et al. 



2001 



2004b 



Zhukov fc Auchere 



20081 ) and whether they are waves at all or rather 



propagating disturbances related to magnetic fie 



with the OME hftoff (e .g. 



Delannee &: Aulanier 



d lin e opening and restructuring associate d 



Wills-Davey et al 



20071 ) 



1999 : 



Ohen et al. 



20021: 



Attrill et al. 



2007 : 



Numerical simulations of EIT waves as fast-mode magnetohydrodynamic al (MHD) 



waves resulted i n wave phenomena rnimicking observatio nal data quite closely (jWang 



2000 



Wu et al. 



2001 



simulations by 



Ofman &: Thompson 



Ofman fc Thompson! (120021 ) and 



2002; 



Ofman 



20071). An intere s ting 



)y-product of the 



Terradas fc OfmanI (120041 ) are stationary 
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brightenings in the vicini ty of the launch sites, in line with observat i ons of such areas at 

J; 



the fronts of EIT wav es (iDelannee fc Aulanier 



Delannee et al 



1999 



Delannee 



200C 



Attrill et al 



20071 ). 



(120071 ) suggested them to be a result of magnetic field restructuring during 



the CME liftoff generated ei ther by Joule 



compression. Simulations by 



Cohen et al 



leating or an increase in density due to plasma 



( 120091 ) favored the latter interpretation with some 



additional effect due to increased temperature resulting from plasma compression. 



Until the launch of the STEREO (Solar- Terrestrial Relations Ob servatory. iKaiser et al. 



20081 ) mission in 2006, with the Extreme Ultra Violet Imager (EUVl; .Howard et al 



20081 ) 



onboard, observations of EIT waves were drastically limited by a ~ 12 min cadence of the 
EIT instrument in the 195 A passband. The identically built instruments EUVI-A and 
EUVI-B onboard the twin STEREO spacecraft observe the solar corona in four different 
EUV passbands with a high observing cadence (up to 75 s) and a large field-of-view (up to 
1.7 R©) from two different vantage points. 

Studies of large-scale waves with STEREO/EUVI reveal velocities in the range 



of ~ 200-400 kms ^ and a decelerating character consisten t with 



large-amplitude MHD fast-mode waves (jVeronig et al. 



Kienreich et al. 



2008 



2010 



reely propagatin g 



Long et al. 



2008 



20091 ). With simultaneous observations by the two STEREO spacecraft, 



it was for t he first time possible t o analyze the 3D nature of EIT waves with stereoscopic 



technique s ( 



Ma et al. 



Kienreich et al. 



2009 



2009 



Temmer et al 



Gopalswamy et al 



Patsourakos fc Vourlidas 



20111 ). 



Veronig et al 



2009 



Patsourakos et al 



(120081), iLong et al 



2009; 



(|2008[ ) and 



(|2009[ ) reported the refraction and reflection of a wave at the border 



of a coro nal hole prov iding strong evidence for the wave nature of the phenomenon. 



However, 



Attrill 



( I2OIOI ) questioned these flndings, favoring instead the hypothesis of two 



simultaneously launched waves from two distinctive l y sep arated sites. In the EUVI wave 



event of 2010 January 17, studied by 



Veronig et al 



(I2OIOI ). flrst observations of the full 3D 
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wave dome and its lateral and radial expansion could be perform ed. In a recent stu dy that 
became available during the revision process of the present paper, 



Long et al 



mm analyze 



the kinematical aspects and the evolution of the full width at half maximum (FWHM) and 
integral of the 2007 May 19 and 2009 February 13 events with similar methods. Therefore 
we will discuss and compare their findings with our results in section |5l For recent reviews 
regarding the kinematics of large-scale EUV waves, their morpholog y an d relationship to 



associ ated solar phenomena we refer to 



fl2010h 



Wills-Davev fc Attrill 



teoioh and 



Gallagher fc Long 



In this paper, we study high-cadence EUV observations of four well defined EUV waves 
observed by STEREO/EUVI that occurred between May 2007 and April 2010. These waves 
are special due to their pronounced amplitudes although they occurred during the extreme 
current solar minimum. Our study has the following main aims: 



1. We apply and compare two different methods to derive the wave kinematics: visual 
tracking of the foremost part of wave front and a method based on the intensity 
profiles of the propagating perturbation. 

2. Based on the perturbation profiles, we derive quantities which provide insight into the 
physical character of the phenomenon in terms of Mach numbers as well as evolution 
of the amplitude, width and integrated intensity of the disturbances. 

3. The perturbation profiles are further studied with respect to associated phenomena, 
such as stationary brightenings, coronal dimmings and rarefaction regions behind the 
wave fronts. 
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2. Data and Observations 

The four large-scale coronal waves under study were recorded by EUVI on 2007 May 19 
2009 February 13, 2010 January 17, and 2010 April 29. The EUVI instrume nts are part of 



the S un Earth Connection Coronal and Heliospheric Investigation (SECCHI; 



Howard et al. 



20081 ) instrument suite on board the STEREO-A (Ahead; ST-A) and STEREO-B (Behind; 
ST-B) spacecraft. The separation angle between ST-A and ST-B steadily increases by 
^45° per year. EUVI observes the chromosphere and low corona in four spectral passbands 
(He II 304 A: T~ 0.07 MK; Fe IX 171 A: T~ 1 MK; Fe XII 195 A: T~ 1.5 MK; Fe XV 
284 A: T~ 2.25 MK ) out to 1.7 Rq with a pixel-limited spatial resolution of 1".6 pixel"^ 



( IWuelser et al. 



20041 ). For the kinematical analysis we used all passbands if possible, while 
for the perturbation profiles only the 195 A images were studied, where the wave signal is 
most prominent. 

The EUVI wave event of 2007 May 19 was recorded by both STEREO spacecrafts 
(ST-A and ST-B) and was observable from the Earth view. High-cadence EUVI images 
in the 171 A (75 s), 195 A (600 sec), 284 A (300 s) and 304 A (300 sec) channels are 



available. This event was the first distinct EUV wave obseryed by 



studied in quite some d e tail by several authors (ILong et al. 



Gopalswamy et al. 



20091 : 



Kerdraon et al. 



2010 : 



Long et al. 



2008; 



]UVI, and was alre ady 



Veronig et al. 



2008 



20111) revealing a propagation 



velocity decelerating from ^ 450 to 200 kms ^ as well as an associated type II burst 
indicating shock formation in the corona. 

The event of 2009 February 13 was observed by both STEREO spacecraft, and 
occurred in perfect quadrature (separation angle of ~ 90°). The wave was observed on disk 
in ST-B and on the limb in ST-A. These observation s proy ided a unique basis to study 



the 3D nature of the wave. 



Patsourakos fc VourlidasI (120091 ) and 



Kienreich et al. 



(120091) 



determined a propagation height of the wave of ~ 100 Mm, with a mean (on-disk) velocity 
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of ^ 230 kms"^. For the analysis, 171 A and 195 A filtergrams of ST-B are available with 
a cadence of 300 s and 600 s, respectively. 

On 2010 January 17, ST-A and ST-B were 134° apart from each other. The wave was 
observed in the eastern hemisphere of ST-B, which was situated 70° behind Earth on its 
orbit around the Sun. This event was unique due to the full 3D wave dome structure that 



was observed in the different EUVI channels with a lateral (on-disk) propagation ve 



ocity 



2010h . 



of ~ 300 kms~^ and an upward propagation velocity of ~ 600 kms~^ (jVeronig et al 

The EUVI-B imaging cadence is 2.5 minutes in the 171 A, 5 minutes in the 195 A, 2.5—5 

minutes in the 284 A, and 5 minutes in the 304 A passband. 

The event of 2010 April 29 was the last and strongest of four homologous waves 
observed by ST- B within a period of 8 hours at a position angle of 70° behind solar limb as 



seen from Earth (JKienreich et al 



1201 ll ). EUVI images in all four wavelengths are available 
during the event. However, the cadence of two hours in the 171 A and 284 A channels 
limited a possible wave tracking to the 195 A and 304 A channels. Due to a lack of wave 
signatures in the 304 A channel, only 195 A filtergrams with a cadence of 300 s are used for 
the analysis. 

All EUVI filtergrams were reduced using the SECCHLPREP routines available within 
SolarSoft. Furthermore, we differentially rotated each data set to a common reference time. 
In order to enhance faint coronal wave signatures, we derived running ratio (RR) images 
dividing each image by a frame taken 10 minutes earlier, as well as base-ratio (BR) images 
dividing each image by the last pre-event image. Finally, a median filter was applied to the 
images to remove small scale variations. RR images are used for the visual tracking of the 
coronal waves, since in RR images the signal of propagating disturbances is highest. The 
perturbation profiles are calculated from BR images, which provide better insight into the 
changes of physical parameters due to the passing wave front. 



3. Analysis 

The kinematical analysis of large-scale propagating disturbances is usually based on 
visual tracking of the outer edge of the wave front. Thus, the results are severely influenced 
by the observer and his/her interpretation of the wave front. Therefore, more objective and 
reproducible as well as automated methods are desirable. Two alternative semi-automated 
methods have been applied so fa r to study EUV waves: the Huygens plotting method 



(IWills- 



2004a 



Davev fc Thompson 



19991) and the perturbation pro f ile method (iWarmuth et al. 



Podladchikova fc Berghmans 



2005 



Muhr et al. 



2010 



Veronig et al. 



2010h 



In this study we focus on the perturbation profile method because of several reasons. 
First, this semi- automated method provides insight into important physical parameters of 
the EUVI wave (amplitude, widths, Mach number). Second, it provides a semi-automated 
alternative to the reconstruction of the wave kinematics by visual tracking. Third, detailed 
analysis of the perturbation profiles provide information on associated phenomena not 
observable by eye, which can give us important insight into the underlying physics of the 
events. In the following, we describe the two different methods, i.e. the visual tracking 
method and perturbation profiles, and their application to the study of EUV waves. 

3.1. Kinematics via visual tracking 



We visually tracked the wave fronts in series of EUVI RR images. The eruption center 
was derived by applying circular fits to the earliest observed EUVI wave fronts o n the 3D 
solar surface appearing as ellipses in the 2D projected image ( j Veronig et al.l 120061 ). In order 
to enhance the statistical significance and to ensure a realistic error estimate for the wave 
center position, all available passbands were used to determine a mean value and standard 
deviation. For each event, we focussed for both, the visual and the perturbation profile 
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method, on the same specific 45° propagation sector, in which the disturbance is most 
pronounced. The chosen sectors are restricted by two great circles (parts of them are shown 
by black curves in Figures [T]-|l]) , which pass through the determined wave center. To obtain 
the wave kinematics, we calculated for each point of the wave front within the selected 
propagation sector its distance from the eruption center along great circles on the solar 
surface, and then averaged over them. 

3.2. Kinematics by perturbation profiles 



The second approach to analyze the EUV wave propagation is based on the perturbation 
profiles. The method starts at the eruption center determined from the visually tracked 
wave fronts (see section 13. ip and sums the intensity values of all pixels between two 
constantly growing concentric circles (again in the deprojected heliospheric plane) defining 
annuli with a radial width of 1° within the selected propagation sector, which span over 
an angular width of 45° (for an illustration see movie 1 of the online version). Due to the 
fact that the area over which we sum up the pixel values is steadily growing by moving 
to greater radii, each annulus is averaged over its pixel sum. This is resulting in a mean 
intensity as a function of distance measured along the solar surface from the wave center. 
The procedure is repeated for each frame, i.e. time step, until the EUV wave fades away. 

In perturbation profiles, propagating disturbances can be identified as a distinct bump 
above the background level (intensity level of 1.0 in base ratio images). Modified Gaussian 
envelopes emphasizing the leading part of the wave bump are then fitted to the perturbation 



profiles for two reasons. First of all, p e rturbation profi 



( jWills-Davev 



2006 : 



Veronig et al. 



2010; 



Kienreich et al. 



es sh ow roughly a Gaussian form 



20111 ) and secondly it is very useful 



for the robustness and automatization of the latter part of t he profile algorithrn , in which 



the positions of the leading edges of the waves are extracted. 



Downs et al. 



fl201l[ ) state that 
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the peak intensity may be strongly affected by the emission of CME material, whereas the 
foremost front corresponds to the actual wave. The trailing edge of the wave is also strongly 
influenced by features like stationary brightenings, flare emission, CME and dimming 
regions behind the wave front. Therefore we use modified Gaussian envelopes to the original 
data set. First of all, the original data is examined for the peak intensity. Then a simple 
Gaussian fit is applied to the first half of the wave bump from the maximum until the 
leading edge. Finally, the Gaussian envelope is mirrored at the position of the maximum. 
The result is a Gaussian fit matching the foremost leading edge of the wave fronts much 
better than simple Gaussian fits while for the rest of the wave bump a characteristic 
envelope shape is formed. 

From the modified Gaussian fits to the perturbation profiles we derive the position of 
the wave front. Starting at the maximum amplitude level of each wave bump we define the 
leading edge of the wave fronts at the position where the Gaussian envelope to the profile 
falls below an intensity level of ///q = 1.02, with I the intensity entries of the current 
image and Iq the intensity entries of the base image. This value is reasonable due to the 
fact that the human eye is able to identify intensity changes of as small as 1-2% above the 
background level. Thus, this should provide us with a good measure when comparing the 
wave front positions from the perturbation profiles and the visual tracking. 

3.3. Calculation of perturbation profile parameters 

It is possible to gain useful information on the wave characteristics by extracting 
distinct wave parameters from the perturbation profiles. We derived the amplitude, the 
FW, the FWHM and the integral below the perturbation profile (down to the intensity 
ratio of 1.02). The FW of the wave is defined as the width of the Gaussian envelope, 
whereas the FWHM is defined as the distance between those two points where the intensity 
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level drops below 50% of the maximuni value. The integral of the perturbation profiles is 
the area below the Gaussian fits. 

The calculation of the magnetosonic Mach numbers is based on the peak perturbation 
amplitude values, Amax, of each wave event, derived from the peak of the Gaussian fits to 
the profiles. Due to the definition of base ratio images, we know that an increase in the 
perturbation amplitude above the background level can be interpreted as A = I/Iq, where 
the intensity / in an optically thin coronal spectral line is given as 



Jh 



n^dh 



'V. 



with the emission measure EM = j^nldh along the line of sight and the c ontribution 
funct ion /(T, rtg) depending on plasma temperature T and electron density rie (jPhillips et al. 
20081 ). To obtain an intensity value in a relative broad EUV passband, like the EUVI 195 A 
passband, the integral over several spectral lines has to be conducted. Since there is a 
temperature and density dependence for each spectral line, the following approximations 
are made: First, th e integration is only along the pressure scale height H with an average 



20061 ). Thus the intensity turns into / ~ f(T,ne)nlH, 



constant density n (iWills-Daveyl 

with nlH the emission measure EM. Secon d, the quiet Sun in th e 195 A passband is 



dominated by its strongest line at 195.12 A ( Del Zanna et al.l 120031 ) which is formed in a 
narrow temperature range of 1.2 MK to 1 . 8 MK , with the peak response at ~ 1.4 MK 



( iFeldman et al. 



19991). In Fig. 3 of 



Zhukovi ( 120111 ) it is evident that the contribution function 



of the 195.12 A spectral line only weakly depends on density, i.e. /(T, Ug) ^ f{T), and only 
weakly depends on temperature near its peak response which is also the position where the 
differential emission measure peaks. Thus and under the assumption that the temperature 
does not change considerably by the passing wave front, the relative intensity change I/Iq 
(with I and Iq referring to the initial a nd final s tate, respectively) can be used to estimate 



the density increase n/uQ ~ (///q)^'^. IZhukovi ( 1201 ll ) gives a detailed discussion of this 
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approach under consideration of the usage of the CHIANTI atomic database. 

For each wave we estimate the density jump X^ = n/riQ at the peak amphtude 
of the perturbation profile. Assuming the shock height to be in the low solar corona 



( IVeronig et al. 



20101 ). we can conclude from typical EUV wave velocities in the range of 



200—500 k ms ^ that the assu mption of a quasi-perpendicular fast-mode MHD shock is 



reasonable ( iMann et al. 



19991 ). Considerin 



perpendicular shock front (e.g., 



Priest 



g the Rankine-Hugoniot jump conditions at a 



19821 ) . the magnetosonic Mach number Mms = Vc/v„ 



(where Vc is the coronal shock velocity and Vms the magnetosonic speed) can be expressed as 

M - / MX. + 5 + 5/3J~ .21 

""' \/(4-^c)(2 + 5/3,/3)' ^'^ 

where Pc is the ambient coronal plasma-to-magnetic pressure ratio, Xc = n/no is the 

density jump at the shock front, and for the specific-heat ratio (the polytropic index) we 

substituted 7 = 5/3. For the plasma beta in th e quiet Sun we assu me a value of (3c = 0.1 



and note that Eq. [2] depends only weakly on (3 (jVrsnak et al 



2002). 



3.4. Coronal dimmings, stationary brightenings and rarefaction regions 

Beside the propagating bright coronal wave front, the perturbation profiles show some 
other distinct features characterized either through persistent intensity depression (coronal 
dimmings and rarefaction regions) or intensity enhancement (stationary brightenings). 

Coronal dimmings are regions of dramatically decreased plasma density occurring after 
the CME lift-off. They are usually interpreted as plasma evacuation due to the field line 



open ing during the CME eruption ( JHudson et al. 



1996 



Harra fc Sterling 



2001 



Harra et al. 



20071 ). and may last for several hours up to a day. We stress that there are different 
types of coronal dimmings, those in the near vicinity of the eruption site, referred to as 
core coronal dimmings, and so-called secondary dimmings further away from the eruption 
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site ( iMandrini et al.l 120071 : iMuhr et al.ll2010l ). In this paper we deal with core coronal 



dimm ings. They usually consist of two areas 



rope (IMandrini et al 



2005 : 



Crooker fc Webb 



decrease of ~ 40-60 % (e.g. 



marking the foot 



2006 



Chertok fc Grechnev 



Muhr et al. 



20051). 



joint s of the erupting flux 



20101 ) and show an intensity 



Stationary brightenings fo rm at the edge of c ore coronal dimmmings and are noticed 



as intensity enhancements (e.g. 



Cohen et al. 



20091 ). They are formed after the wave passage 



and last for a period of several tens of minutes until fading away. A possible explanation for 
their formation is that they are a result of material flow, storage an d compression at this 



speci f ic location due to an obstacle that can not be overcome easily ( iDelannee fc Aulanier 



1999 



Delannee et al. 



20071 1. 



Rarefaction regions are expected to occur right behind a propagati ng wave pulse, 



clearly visible by an intensity amplitude below the background intensity (JLandau &: Lifshitz 



19871 ). The moving wave pulse is a region of compressed plasma, continuously compressing 
the plasma in front of itself. Consequently, in the section behind the wave pulse a rarefaction 
region is formed, expanding to a certain distance. The wave travels at a velocity higher than 
the surrounding Alfven velocity and thus the region behind the pulse is thinned out. The 
pressure and density in this region falls below the equilibrium values. The rarefaction region 
is the result of this evacuation due to the wave propagation. In contrast to the coronal 
dimming, which may be stationary and in most cases relatively huge, the rarefaction region 
propagates following the wave pulse and can be noticed as a small, localized dip in the rear 
section of the perturbation profile between the wave pulse and the coronal dimming. 
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4. Results 

4.1. Wave kinematics 

In Figures [IHll the morphology and evolution of the four EUVI wave events under study 
is shown in BR images. We note that the wave fronts are not observed in the near vicinity 
of the eruption center but are first visible at distances of ~ 240 Mm, 90 Mm, 150 Mm and 
260 Mm from the respective initiation centers of 2007 May 19, 2009 February 13, 2010 
January 17, 2010 April 29 event and can be followed up to distances of 600-900 Mm. 

The top panels of Figures [SH5] show the distance-time plots of the events. We display 
the wave kinematics derived from both methods (visually tracking and perturbation 
profiles) together with error bars, linear and quadratic least square fits to the visually 
tracked wavefronts, as well as 95% confidence intervals and prediction boundaries of the 
linear fits. The mean velocities derived for the four events under study lie in the range of 
220-350 km s-^ 

Three of the four events propagate with constant speed. The small deceleration values 
derived from the quadratic fits are in the range of —6 to — 13 m s~^, which corresponds to a 
deceleration of ~ 5% with respect to the start velocity values derived from the quadratic fit 
to the wave kinematics. This lies within the confidence intervals, and thus it is justified to 
represent the kinematics of the three events by a propagation with constant velocity over 
the full propagation distance. We note that the kinematical analysis of the perturbation 
profiles of these three events reveal a somewhat stronger deceleration as compared to 
the visual tracking method, but they are still only marginally significant corresponding 
to a deceleration rate of ~ 10% of the start velocity values. Only for the 2007 May 19 
event, we obtain a significant deceleration of —193 m s~^, which corresponds to 30% of 
the start velocity. This is also the fastest wave in our sample with a mean velocity of 
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f« 350 kins"^. In Table l^?T| the derived information on the propagation characteristics 
(velocity, deceleration) are summarized. 

Event S/C vun [kms"^] Vquad.o [kms^^] ao [m S"^] /max Xc Afms 

vis. prof. vis. prof. vis. prof. 

2007 May 19 ST-A 348±29 338±18 483±73 429±58 -193.2 ±87 -84.7 ± 60 1.61 1.27 1.21 

2009 Feb 13 ST-B 228±17 221±36 237±72 261±72 -6.9 ± 25 -23.4 ± 30 1.39 1.18 1.14 

2010 Jan 17 ST-B 286±7 275±17 300±24 314±53 -13.4 ±38 -39.4 ±49 1.59 1.26 1.20 
2010 Apr 29 ST-B 337±27 321±23 342±88 361±76 -7.6 ±52 -51.7 ±93 1.22 1.12 1.08 

Table 1: EUVI wave properties of the events under study. We list the mean velocity values 
viin derived from linear fits, as well as the start velocity values Vquad,o and the acceleration 
values ao derived from the quadratic fits for both methods, i.e. visual tracking method and 
profile method. Additionally, the peak intensity-amplitudes extracted from the perturbation 
profiles together with the density jump values Xc and the calculated Mach numbers Mms are 
listed. 

We note that the position measurements of both methods match each other quite well, 
the maximum differences lie within an average of 20-40 Mm. As a consequence, also the 
velocity values derived via visual tracking correspond to these obtained by the perturbation 
profiles, and both values agree within the error limits (see Table S]T]). 



4.2. Overall characteristics of perturbation profiles 

The consistent kinematical results obtained with both methods suggest the perturbation 
profile method to be an adequate alternative to the visual tracking method. In particular, 
it is an additional possibility to analyze large-scale disturbances in the solar atmosphere 
by producing kinematical results not burdened by subjective judgements of the observer. 
Thus, the perturbation-profile results are easier to reproduce than those derived by the 
visual-tracking method. 
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Moreover, we can quantify various important wave parameters from the perturbation 
profiles. Figures [n]-[I2] show for the four wave events under study the evolution of the 
perturbation profiles. In all four cases we observe a clear intensification in their early phase 
reaching its peak amplitude Amax value around 10-20 minutes (second or third panel of 
Figures [9]- [T2l) after the first remarkable wave bump can be observed. The evolution of 
the intensity amplitudes of all events under study is displayed in the middle panels of 
Figures ISHHl The peak perturbation amplitudes Amax lie in the range of 1.2 — 1.6. The 
calculated density jump values, Xc, are therefore in the range of 1.1—1.3, and the peak 
magnetosonic Mach numbers, M^s, lie in the range of 1.08—1.21 (Table H7I]) . 



In addition, we derived the full width of the wave as well as the FWHM. The evolution 
of the pulse width is plotted in the bottom panel of Figures [5] - [8] for the four wave events 
under study. The width of the wave pulse as well as the FWHM increases during its 
evolution by a factor of 2-3, whereas the integral below the perturbation profile remains 
basically constant for the events of 2007 May 19 and 2010 April 29, and decreases for the 
events of 2009 February 13 and 2010 January 17 by a factor of 4 and 3, respectively. 

4.3. Distinct features observable in the perturbation profiles after the wave 

passage 

The inspection of the derived BR perturbation profiles reveals a diversity of distinct 
features after the wave passage, including coronal dimmings, stationary brightenings, and 
rarefaction regions. In the following, these are discussed individually for each of the four 
events in our sample. 
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4.3.1. Event of 2007 May 19 



The profiles of 2007 May 19 sliow tliree distinct pliases: tlie leading wave pulse 
(enveloped by gaussian fits), a reversal point (where the intensity level switches from values 
greater than 1.0 to below 1.0) and the trailing rarefaction region. Additionally, in the rear 
section of the perturbation profiles we observe a deep core dimming extending from the 
initiation site up to ~ 300 Mm, with an intensity decrease of ~ 70% lasting for several 
hours (Fig. [H]). 

The rarefaction region and its propagation is of special interest. This phenomenon is 
theoretically predicted to be located behind the propaga ting front of a compression wave 



but observations are barely reported (I Cohen et al 



20091 ). From 12:52 UT until 13:02 UT 



we observe the formation of a deep coronal dimming behind the wave pulse. During its 
early evolution it becomes more pronounced and reaches its maximal spatial extension at 
13:02 UT. After 13:02 UT it is quasi stationary. However in the frontal part of the dimming, 
from 300 Mm up to the trailing part of the wave pulse, a smaller, propagating dip is clearly 
visible. We interpret this as the rarefaction region following the wave front. Its formation 
is a combination of two aspects: 1) the driver stops at a distance of ~ 300 Mm (where the 
outermost point of the deep-dimming region is located), and 2) behind the wave pulse the 
wave forms a region where the density and pressure fall below the equilibrium level. 

This event was associated with a metric type II burst, indicating shock formation in 
the corona. We use the observations of the type II burst to obtain an alternative estimate 
of the shock Mach number. Type II bursts usually show the fundamental and harmonic 
emission band, both fre quently being split in two parallel lanes, so-called band-split 



(INelson fc Melrose 



19851 ) . The interpretation of the band spli t as the plasrna emis sion from 



the upstream and downstream shock regions was affirmed by 



Vrsnak et al. 



mm . Thus, 



the band-split can be used to obtain an estimation on the density jump at the shock front 
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and thus the Mach number of the ambient plasma (jVrsnak et al.ll2002l ) 



The composite dynamic radio spectrum (Fig. [T5]) shows complex and intense radio 
emission, consisting of a group of type III bursts and a type II burst. For a detailed 
dy namic spectrum analy sis covering the frequency range 0.4-300 MHz, we refer to Fig. 4 



m 



Kerdraon et al. 



( I2OIOI ). In the present study we focus on the type II burst that started 
around 12:51:30 UT at a frequency of 160 MHz. In the period 12:51:30 UT to 12:54:00 UT 
we were able to recognize a band-split pattern (indicated in Fig. [13] as two black lines), with 
a relative bandwidth of BDW = A/// ^ 0.19-0.22. The relative bandwidth BDW is 
determined by the density jump at the coronal shock front Xc = P2c/pic, where_£L^^_and_£2 



are t he densities upstream and downstream of the shock front (for details see 



Vrsnak et al. 



200l|). Since BDW= (/a - /i)//i = aZI^^^iJ - 1. we find X, = 1.41 - 1.48. Using a 



5-fold Saito coronal density model (jSaito 



19701), we find from the derived values of the 



density jump Xc a Mach number of 1.32-1.40. 

From the peak amplitude of the EUVI perturbation profile we derived Mms ~ 1.21. A 
possible explanation for the smaller estimate is that, due to the width of 45° over which the 
perturbation profiles were derived, eventual higher values are averaged out (see Figure [9]). 
The curve displayed by small black crosses shows the overall 45° sector we used for our 
analysis, while the grey line profile is displaying a small 2° sub-sector in the central part of 
the overall 45° sector. By shrinking the sector to a narrower width around the most intense 
propagation direction, we obtain a considerably higher peak due to a lower signal averaging, 
and the Mach number is increased to values of at least Mms = 1-3. We note that the sector 
narrowing and the subsequent increased Mach number implies that the extracted values 
are (in fact always) underestimated and thus a lower limit for the Mach number estimate. 
Comparing the results of both approaches we find the Mach number of the ar nbient plasma 
to be Mms = 1.2-1.4, a reasonable range of values for a coronal shock wave ( iVrsnak et al 
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200ll ). Since EUV waves and type II bursts are generated at different heights in the solar 
atmosphere the Mach numbers are not necessarily the same. In particular, type II bursts 
are point-like sources formed where the disturbance is most intense, corresponding to higher 
Mach numbers. Shocks having Mach numbers in the range of 1.2-1.4 as determined for the 
events under study are usually considered as weak shocks. 



4.3.2. Event of 2010 January 1 7 



The perturbation profiles of 2010 January 17 show persistent stationary brightenings 



(Figures [Toll . They appear after the waye fron 



( iDelannee fc Aulanier 



1999 



Attrill et al 



2007 : 



; passage and are on ly slowly fading away 



Delannee et al 



20071). The EUV wave starts 



at 03:56 UT as a relatively small wave pulse that is intensifying until 04:01 UT. In the 
subsequent panels of Figure [TOl after the wave has intensified up to its maximum amplitude 
of 1.6 at 04:01 UT, a strong and long-lasting, quasi-stationary bright feature becomes 
clearly visible. We find first evidence for this stationary parts at 04:06 UT at a distance of 
100-300 Mm from the source location, i.e. at the distance of the first wave front appearance 
(see Figure fTOl) . It is well pronounced and exceeds an amplitude level of 1.25, lasting for 
at least 30 min after the wave passage. Compared to the time scale of the propagating 
disturbance it is only slowly fading away. This evolution can also be seen in the BR images 
in Figure [3l where we overplotted the position of the stationary brightenings by grey lines. 



4.3.3. Event of 2009 February 13 



The perturbation profiles of the quadrature event of 2009 February 13 are shown in 
Figure [Til The evolution of the most enhanced direction of 45±22.5° is displayed in the 
perturbation profiles by crosses. The propagating wave pulse can be clearly identified. 
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The grey line profile represents a different sector with its main direction to be 270±22.5°. 
It is characterized by an intensity enhancement which can be interpreted as a stationary 
brightening, showing a similar evolution like the one of 2010 January 17. The maximum 
intensity of these stationary brightenings is even more pronounced than that of the 2010 
January 17 event and exceeds a level of 1.35 at 05:55 UT at 100-300 Mm from the source 
location (see Figure [TTl) . It is again induced by the wave front passage and only slowly 
fading away. The propagating wave pulse ahead of the stationary brightening is not as 
pronounced as for the main direction but still detectable although it is considerably smaller 
than the stationary brightening. In both directions the coronal dimming is clearly visible. 
Its location is restricted to the area behind the stationary brightening at a distance up to 
^ 100 Mm. 

4.3.4. Event of 2010 April 29 
The wave event of 2010 April 29 is the f ourth and strongest o f four homologous waves 



that occurred within a time range of 8 hours (JKienreich et al. 



2OIII ). The profile evolution of 



the event is displayed in Figure O We can identify three different parts in the rear section 
of this perturbation profile: a dimming region, a stationary brightening and a rarefaction 
region, all visible after the waves' passage. The dimming region is prominent from the first 
image at 06:20 UT until the last profile shown at 06:40 UT, occurring at 0-100 Mm from 
the source location. Right in front of it an intensity enhancement is present. A definite 
identification with a stationary brightening is difficult due to its relatively weak appearance 
compared to the background level of 1.0. This brightening lasts for at least 40 min after the 
passage and is only slowly fading away. In front of this intensity enhancement a dip in the 
profile is propagating away from the initiation center, which we interpret as the rarefaction 
region, similar to the event of 2007 May 19. 
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5. Discussion and Conclusions 

1. We analyzed four well pronounced EUV wave events observed by the STEREO 
EUVl telescopes in order to compare two different kinematical analysis techniques, 
the generally used visual tracking method and the semi-automated perturbation 
profiles. The differences in the determined positions of the leading edge of the 
wave fronts using both methods are maximal 40 Mm. Thus, we conclude that the 
perturbation profiles are a suitable method to analyze large-scale waves. The big 
advantage of the perturbation profile method is the higher degree of automatization 
and reproductiveness due to the objective measurements of the wave location. 
Nevertheless, there are some restrictions. In the later evolution phase, the wave fronts 
become more irregular showing changes in shapes and propagation direction. These 
effects are not considered in the profile method. Thus, due to the lower amplitude 
and the smearing of irregular fronts the wave profile is no longer distinct against the 
background. These properties lead to a systematic underestimation of the distance 
of the wave compared to the visual tracking method. Additionally, the diffuse bright 
fronts are easier to identify in RR images due to the fact that they have a higher 
contrast than BR images. 

2. The determined propagation velocity derived from both methods match within 
the error limits. For the four wave events under study, we obtain mean velocity 
values in the range of 220 — 350 kms^^, i.e. within the velocity rang e for fast 



magnetosonic waves during quiet Sun conditions (e.g. 



Mann et al. 



19991). Three of the 

four kinematical curves show only small deceleration values of —6 up to — 13 m s^^, 
which lie within the error range. For the fastest event of 2007 May 19, the obtained 
deceleration of —193 m s~^ is significant. 



Long et al 



(l201l[ l has analyzed the kinematics of the May 19, 2007 and February 
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13, 2009 events with a similar perturbation profile method. They find for the May 
19, 2007 event the following values: start velocity fo=447±87 km s~^, acceleration 
a=-256±134 m s~^, a broadening of the pulse width from 50 to 200 Mm and a 
decrease of the peak amplitude intensity from 60 to 15% of the background values. 
Our findings are fo=429±58 km s~^, a=-85±60 m s~^, and similar results for the 



broadeni ng of the pu 



13, 2009, 



Long et al. 



se ana 



their intensity evolution. For the second event of February 
fl201lh find t;o=274±53 km s-\ a=-49±34 m s'^ while our 
results are fo=261±72 km s~^, a=-24±30 m s~^. The profile broadening and the 
intensity evolution is again similar. For the February 13, 2009 event the outcomes of 
both studies are consistent within the error limits. The results of May 19, 2007 show 
differences in the acceler ation values which can be explained by taking into account 



that the method used by 



Long et al. 



(J201ll ) and ours is not identical. 1) While they 



are extracting the position of maximum intensity we are using the leading edge of the 
wave fronts for our analysis. The difference in the acceleration values is thus due to 
this usage of different features for the position measurements. The maximum intensity 
peaks are propagating at a slower speed compared to the leading edges, an effect that 
is expected for a feature that is broadening during its evolution. Due to that the 
deceleration values are smaller for the leading edge measurements while the st a rting 



Long et al 



(1201 ll ) 



velocity values are not affected. 2) The overall sector width used in 
is varying from event to event while we are using a constant angular width of 45° in 
each event. 3) The main direction of the calculated sector as well as the the initiation 
centers are not exactly the same (the uncertainty is about 20 Mm). 

3. Considering the events under study as low-amplitude MHD fast-mode waves the 
Mach numbers and wave velocity values are proportional to each other, Mms = v/vms- 
Thus, from the derived Mach number evolution (from the maximum value down to 



Mms ~ 1) we expect a decrease of the propagation velocity by Av ^ 100 kms 



-1 
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(2007 May 19), ^ 35 kms"^ (2009 February 13), ^ 60 kms^^ (2010 January 17) and 
p^ 30 kms~^ (2010 April 29). The velocity changes we derive from the least-square 
quadratic fits reveal ^ 200 kms^^ (2007 May 19), ^ 40 kms^^ (2009 February 13), 
^ 20 kms-i (2010 January 17) and ^ 30 kms'^ (2010 April 29). For the 2009 
February 13 and 2010 April 29 events, the derived velocity changes are in the same 
order of ~ 30 kms~^ as the error on the velocity determination. Hence, the weak 
deceleration is hidden in the measurement uncertainties. For the 2007 May 19 event 
the result of 200 kms~^ clearly exceeds the error in velocity of ^ 70 kms^^, thus the 
deceleration is observable and significant. According to the Mach number evolution 
for the 2010 January 17 event, an observable deceleration of the propagating wave 
is expected. The determined velocity change of ~ 60 kms~^ exceeds the velocity 
error of 20 kms~^, thus masking of the deceleration by data scatter is unlikely. 
Nevertheless, from the kinematical measurements obtained with the visual tracking 
and the perturbation profiles we do not derive a significant deceleration, as it would 
be expected from the peak Mach number estimated for this event. 



4. Each event can be followed over a period of at least 30 minutes, during which we first 
observe an intensification of the wave, which is followed by a steady decrease and 
broadening of the wave pulses as extracted from the perturbation profiles. From this 
we derived useful information on the characteristics of the wave pulses, like the peak 
amplitude, the width, and the integral of the wave pulse. The evolution of the pulse 
width shows a clear broadening over time by a factor of 2-3, similar t o the broadening 



of the 



2004bl : 



WHM values. These values are typical for large-scale waves (jWarmuth et al 



Wills-Davev 



2006 



Warmuth 



2010 : 



Veronig et al 



20101 ). The integral below 



the perturbation profile basically remains constant (2007 May 19 and 2010 April 29) 
or decreases by a factor of 3-4 (2009 February 13 and 2010 January 17). We want 
to stress that the combination of profile broadening and amplitude decay leading to a 
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constant integral below the disturbance profile is co nsistent with the characteristics of 



a freely propagating wave (JLandau fc Lifshitz 



19871 ) 



5. The peak magnetosonic Mach numbers derived for the events under study are in 
the range of Mms = 1.08-1.21, indicative of the evolution of weak coronal shocks. 
The magnetosonic Mach number for 2007 May 19 was determined via two different 
approaches leading to similar results in a range of M^s = 1.2-1.4. Thus, the usage of 
the perturbation maximu m amplitude for the dete rmination of the magnetosonic Mach 



number seems reasonable. 



Grechnev et al. 



(l201l[ ) discussed the problem of projection 



effects in measurements of EUV wave. Since we do not observe the perturbation from 
a view-point that is located directly above the propagating disturbance but from an 
inclined position, the intensity values are always underestimated. This is probably 
also the reason why we do observe Gaussian wave pulse profiles instead of a sharp 
edge, expected in the case of a shocked disturbance. 

6. We observe three different distinct features in the perturbation profiles occurring 
after the wave passage. In three of four events deep coronal dimmings are prominent 
at small distances from the eruption center. Only for the 2010, January 17 event 
it is hardly present in the perturbation profiles of the used sector. However, an 
inspection of Figure [3] reveals the existence of coronal dimmings in the near vicinity 
of the initiation site most prominent to the western sectors close to the limb and also 
off- limb. 

The second distinct feature in the wake of three of the four the EUV waves are 
stationary brightenings. These stationary parts appear after the waves' passage at 
distances of 100-300 Mm and fade only slowly away. Due to the fact that the the most 
intense wave profile always forms one frame before the stationary brightening appears, 
these parts seem to be an obstacle for the wave front. These stationary brightenings 
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are not circumferential around the eruption center but confi ned to specific sect ors 



CohenetaL 



(120091) 



which is best observed for the event of 2009 February 13. 

did numerical simulations on this event and showed that stationary brightenings 

are developing at the outer edges of the core coro nal dimmings. Obsery ations by 



Delannee &: Aulanierl f|l999f 



simulations by IChen et al. 



De 



anned ( 120001 ) and 



Delannee et al. 



(120081) as well as 



fcOOSl l indicate that these stationary brightenings appear 
at the locations where the connectivity of the magnetic field lines changes. Since 
these stationary brightenings are located in front of the deep core dimming regions 
(interpreted as the footprints of the expanding CME), are restricted in their expansion 
(their maximal distance is given by the first wave fronts location) and appear only in 
confined sectors after the waves' passage, they can be interpreted as a signature of the 
CME expanding fianks. 

The third distinct feature are rarefaction regions observable in two of four wave 
events, on 2007 May 19 and 2010 April 29. The rarefaction region behind the wave 
pulse is a prospective fea ture that develops trailin g a large-amp li tude perturbation 



(jLandau fc Lifshitz 



19871 ). Indeed, simulations by 



Cohen et al 



f l2009h revealed 



regions, which are dim in intensity and propagate after the bright wave front. 



The EUV wave velocities are in the range of fast magnetosonic waves in the quiet 
solar corona and show either a decelerating characteristics or constant velocity (within 
the measurement uncertainties) during propagation. The perturbation profiles reveal 
intensification, broadening and decay of the propagating wave pulse, while the integral 
evolution is either constant, or decreasing. The small Mach numbers are indicative of 
a weak coronal shock. Stationary brightenings formed in confined sectors at the outer 
edge of the core coronal dimming can be interpreted as signatures of the expanding CME 
flanks. The first observable wave fronts are always located right in front of them. Thus, we 
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conclude that the EUV wave events can be interpreted as freely propagating, weak shock 
fast-mode MHD waves initiated by the CME lateral expansion. 
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Fig. 1.— Sequence of STEREO/EUVI-A BR images of the 2007 May 19 event. The deter- 
mined initiation center for the EUVI wave is indicated by the cross. The black hnes indicate 
great circles through the wave center determining the propagation direction of 10° ± 22.5°, 
on which we focus the wave analysis. The black fronts correspond to the positions of the 
leading edges of the wave fronts extracted from the perturbation profiles. The field of view 
(FoV) is x=[-1200", +1200"], y=[-1200", +1200"] with the origin at the center of the Sun. 
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Fig. 2. — Same as Figure[T]but for the 2009 February 13 event, with the propagation direction 
of 340° ± 22.5°. 
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Fig. 3. — Same as Figure[3]but for the 2010 January 17 event, with the propagation direction 

of 65° ± 22.5°. 
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Fig. 4. — Same as Figure H] but for the 2010 April 29 event, with the propagation direction 
of 300° ± 22.5°. 
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Fig. 5. — Top: Kinematics of the wave fronts of the 2007 May 19 wave observed on the solar 
disk in the EUVI-A 17lA (triangles) and 195A (diamonds) channels by the visual tracking 
method together with the kinematical measurements we derived from the perturbation pro- 
files (full circles). The error bars added to the measurements represent the diffusiveness of 
the wave fronts. The solid and dashed lines indicate the linear and quadratic least square 
fits to the data set obtained by visual tracking, respectively. The dotted lines indicate the 
95% confidence interval (filled grey area) and the prediction boundary of the linear fit. Mid- 
dle: Evolution of the amplitude of the perturbation profiles shown in Figure O Bottom: 
Evolution of the frontal width of the wave pulse (triangles), the full width at half maximum 
(FWHM, squares) and the integral of the gaussian envelopes (asterisks). 
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Fig. 6. — Same as in Fig. |5]but for the 2010 January 17 wave. Top: Kinematics of the wave 
fronts are derived from the EUVI-B 171A (triangles), 195A (diamonds), 284A (squares) and 
304A (asterisks) observations. 
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Fig. 7. — Same as in Fig. [5] but for the 2009 February 13 wave. 
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Fig. 8. — Same as in Fig. [5] but for the 2010 April 28 wave. 
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Fig. 9. — Perturbation profiles of 2007 May 19 ST-A observations in tlie analyzed sector 
(10° ± 22.5°) togettier witli tlie gaussian envelopes of the wave pulses (thick grey curves). 
The thin grey perturbation profiles correspond to the narrow sector of 8±2°. The grey arrows 
point to the core dimming while the black arrows point to the rarefaction region. 



-36- 



03:51:02 UT 



i^^rt ll ^^ ll ^W l ^ llllll h 



03:56:02 UT 




J j^^^^ lllll W ll H^ I M l W I 



04:01:02 UT 




04:06:02 UT 




04:1 1:02 UT 



+ 




i»^. 



^ I ^ M^^ 4 



04:13:32 UT 



^^^*+^ 



^ 



^4^^ 



200 400 600 

Distance [Mm] 



Fig. 10. — Perturbation profiles of 2010 January 17 ST-B observations in tlie analyzed sector 
(65° ± 22.5°) together with the gaussian envelopes of the wave pulses (thick grey curves). 
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Fig. 11. — Perturbation profiles of 2009 February 13 ST- A observations in tlie analyzed 
sector (45° ± 22.5°) together with the gaussian envelopes of the wave pulses (thick grey 
curves). The thin grey perturbation profiles correspond to the sector of 270° ± 22.5°. 
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Fig. 12. — Perturbation profiles of 2010 April 29 ST-B observations in the analyzed sector 
(65° ± 22.5°) together with the gaussian envelopes of the wave pulses (thick grey curves). 
The black arrows point to the rarefaction region. 
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Fig. 13. — A composite dynamic r adio spectrum in the frequency range 0.4-300 MHz 



(adapted from 



Kerdraon et al 



(I2OIOI )). Type II and type III bursts are clearly visible. The 



two branches of the bandsplit of the fundamental emission band are marked by black lanes. 
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